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By using a newly developed di-nuclear system model with a dynamical potential energy surface— 
the DNS-DyPES model, hot fusion reactions for synthesizing superheavy nuclei (SHN) with the 
charge number Z = 112-120 are studied. The calculated evaporation residue cross sections are in 
good agreement with available data. In the reaction °°Ti+?“9Bk — 7°°-"119 + an, the maximal 
evaporation residue (ER) cross section is found to be about 0.11 pb for the 4n-emission channel. 
For projectile-target combinations producing SHN with Z = 120, the ER cross section increases 
with the mass asymmetry in the incident channel increasing. The maximal ER cross sections for 
*8Fet+?4Pu and °4Cr+4*Cm are relatively small (less than 0.01 pb) and those for °°Ti+?4°Cf and 
50Ti+?51Cf are about 0.05 and 0.25 pb, respectively. 


PACS numbers: 24.10.-i, 25.70.Jj, 24.60.Dr, 27.90.+b 


In the last decades, a lot of experimental progresses 
have been made in synthesizing superheavy elements 
(SHE). Up to now, SHEs with charge number Z < 118 
have been produced via cold fusion reactions with Pb or 
Bi as targets [1, 2] and hot fusion reactions with *%Ca as 
projectiles [3, 4]. There have been also some attempts 
to synthesize superheavy nuclei (SHN) with Z > 118. 
For example, experiments with projectile-target combi- 
nations °°Fe+?44Pu [5] and °°Ti+?4°Cf [6] have been 
performed in order to produce the element 120 but no 
a decay chains consistent with fusion-evaporation reac- 
tion products were observed. 

The evaporation residue (ER) cross section opr of fu- 
sion reactions depends strongly on the projectile-target 
combination and the incident energy. The study of such 
dependences is interesting and useful particularly when 
one tries to synthesize new SHEs with Z > 118 because 
Ogr Of reactions with these nuclei as evaporation residues 
becomes tiny which makes the experiment much more 
difficult. In recent years much effort has devoted to the 
investigation of the synthesis mechanism of SHN with 
Z > 118. Using a dinuclear system (DNS) model, Feng 
et al. calculated the cross sections of cold fusion reac- 
tions with isotopes of elements 119 and 120 as evapora- 
tion residues and the maximal opr was predicted to be 
about 0.03 pb and 0.09 pb for elements 119 and 120, re- 
spectively [7]. Hot fusion reactions with °°Ti as the pro- 
jectile were studied extensively by using the DNS mod- 
els [8-11], the fusion by diffusion models [12-15], and 
some other models [16, 17]. The optimal incident energy 
and the maximal opr from different models or with dif- 
ferent parameters vary very much. For example, cgr for 
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the reaction °°Ti+74°Cf > 296120 + 3n ranges from 1.5 
to 760 fb [11]. 

In this work, we will study hot fusion reactions pro- 
ducing SHN with Z > 112 and in particular the elements 
119 and 120 using a newly developed DNS model with a 
dynamical potential energy surface (DyPES) [18, 19] (the 
model is termed as the DNS-DyPES model). The impor- 
tance of dynamical deformations of fragments in dissi- 
pative heavy-ion collisions has been long known [20, 21]. 
The dynamical deformations of the projectile and the 
target in the entrance channel have been included in the 
DNS models [22-24] and they are crucial for calculating 
the local excitation energy of a DNS during the process 
of nucleon transfers. Quite recently, Huang et al. [25] 
developed a new DNS model which takes into account 
the dynamical deformations of each DNS. In that model, 
a three-dimensional master equation is solved with three 
variables, the deformations 8; (i = 1 and 2 for each nu- 
cleus in a DNS) and the mass asymmetry 7. In order 
to i) take into account the influence of the dynamical 
deformations of nuclei in each DNS and ii) make the cal- 
culation easier, we treat the dynamical deformations in a 
more transparent and economic way in the DNS-DyPES 
model. The details of this model will be published else- 
where [19] and here we only discuss briefly how to include 
the dynamical deformation. 

As usual, the evaporation residue cross section in a 
heavy-ion fusion reaction is calculated as the summation 
over all partial wave J, 


opR(Ec.m.) = 5 Ocap (Beina J) Pen (Eem.; J) 
J 


x Weur(Ec.m., J), (1) 


with E.m, the incident energy in the center of mass 
frame. In this work the capture cross section Ocap 
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is calculated with an empirical coupled channel ap- 
proach [26, 27] and the survival probability Wsur is calcu- 
lated using a statistic model [28, 29]. In the fusion pro- 
cess, an excited compound nucleus (CN) may be formed. 
We calculate the formation probability of a CN Pon 
based on the DNS concept. The DNS concept was pro- 
posed by Volkov in order to describe the deep inelas- 
tic transfer process [30]. This concept was later used 
to study the competition between complete fusion and 
quasi-fission and to calculate the fusion probability in fu- 
sion reactions [31-33]. Based on the DNS concept, sev- 
eral models have been developed for the study of the 


Pet 


= Wa, ar (t 


Since Ay+A2 = Ap+Ar, only A; is explicitly included in 
the above equation. d4, (t) is the microscopic dimension 
for a DNS (Aj, A2) with a local excitation energy E ys 
defined in Eq. (3). Ejjng is shared by the two nuclei in 
this DNS. For each nucleus, a valence space is opened due 
to the excitation and those nucleons in the states within 
the valence space are active for the transfer between the 
two nuclei. da, (t) = Cmi CN2 where Ny is the number 
of valence states and mz, is that of valence nucleons [22]. 
A is the quasifission probability of the DNS (Aj, A2) 
and Wa, a; (t) = Wa, a, (t) is the mean transition proba- 
bility between the DNS’s (Aj, Az) and (A1, A$). For the 
details about how to solve the master equation, please 
refer to Refs. [22-24]. Here we only focus on the local 
excitation energy and the dynamical deformations. 

In Eq. (2), da, (t), A% , and Wa, 4; (t) all depend on 
the local excitation energy of the DNS, 


Epns(A1,t) = Erota — Epyns (41, t) — Epns(t), (3) 


with 
Etotal = Eom. =+ (Mr + Mp)c’, (4) 
Epns(A1,t) = Vons(A1,t) + (Mi + M2), (5) 
i J(J+1 
Epns(t) = ( ) (6) 


2Jpns(41,t) 


Vpns(A1, t) = Vn(Ai,t) + Vo(A1,¢) and Vn(Aı,t) and 
Vo(A1,t) are the nuclear and the Coulomb interactions 
between the two nuclei. The potential energy in the mass 
asymmetry degree of freedom, which is often called as the 
driving potential at t = 0, is defined as 

Vers(A1, t) = Vons(A1, t) + (Mı + Mə — Mr — Mp)c’. 
(7) 
The interaction potential Vpns(A1,R,t) between the 
two nuclei in a DNS depends on the distance be- 
tween their centers R and Vpns(Ai,t) in Eqs. (5) 
and (7) takes the minimum value of the pocket with 


) [da (t)P(AY, t) — da; (E) P(A) 


synthesis mechanism of SHN (see [34-37] and references 
therein). 

The basic idea of the DNS concept is that after the 
capture process, a DNS (Ap, Ar) in the entrance chan- 
nel is formed. Then the DNS evolves via nucleon trans- 
fer along the mass asymmetry coordinate 7 instead of in 
the direction of the relative distance between the pro- 
jectile and the target R. During the nucleon transfer 
process, any DNS (Aj, Ag) with A; = 0,1,---,Ap + Ar 
and Ag = Ap + Ar — A; may be formed. The evolution 
of the distribution function of each DNS with time can 
be described by a master equation [22-24], 


— AY (t)P(A1,)- (2) 
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FIG. 1. (Color online) The total “intrinsic” energy Eint (the 
black solid curve) and the potential energy Vprs (the red dot- 
ted curve) as a function of the dynamical deformation 66 for 
the projectile-target combination °°Ti + 7“°Cf. The vertical 
dashed line shows the maximal dynamical deformation Bmax. 
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respect to R in Vpns(A1, R,t), 
Vpns (Ai, R, t)|R=Rpocket £ 

Due to the attractive nuclear force and the repul- 
sive Coulomb force, both nuclei in a DNS are distorted 
and dynamical deformations develop during the process 
of nuclear transfers [20, 21]. This results in the time- 
dependence of the potential energy surface (PES). The 
nuclear interaction is calculated with a double-folding 
method [38] and the Coulomb interaction from the Wong 
formula [39]. In this work we assume a tip-tip orienta- 
tion of the two deformed nuclei and that the dynami- 
cal deformations of the two nuclei satisfy 5820, /A, = 
583C2/Az2 [26] with the stiffness parameter C; (i = 1 and 
2) calculated from a liquid drop model [40]. We then 
define 68 = (68, + 662)/2 and following Refs. [20, 21] 
we assume that the dynamical deformation evolves in an 


1.€., Vpns(A1, t) = 
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FIG. 2. (Color online) (a) The maximal dynamical deforma- 
tion 6Bmax and (b) the dynamical potential energy surface 
(DyPES) defined in Eq. (7) at t = 0 and t = oo as functions 
of the mass asymmetry coordinate 7 for the projectile-target 
combination °°Ti + 74°Cf. The vertical dashed line shows the 
entrance channel. 


overdamped motion, 
dB(t) = 68 max (1 = galaa) , (8) 


where the relaxation time Taget = 40 x 10~?? s [21] and 
the maximal dynamical deformation Bmax is determined 
by minimizing the total “intrinsic” energy, 


Fint(A1, 68) = Vn(A1; 81, B2) + Vo(A1; b1, 62) 
i 


i=1,2 


where the quadrupole deformation 6; = 6? + 68; (i = 1 
and 2) with the static deformation parameters 8? taken 
from Ref. [41]. This is illustrated for °Ti + 74°Cf in 
Fig. 1 (the black solid curve) where one finds that with 
68 increasing, Fin, decreases and takes the minimal value 
at 66 ~ 0.36. 

Figure 2 shows the maximal dynamical deformation 
OPmax and the DyPES defined in Eq. (7) at t = 0 and 
t = oo as functions of the mass asymmetry coordinate 
T = (Ay = Apo) /(At + A2) for Opi + 249 Of, The val- 
ues of nuclear masses are taken from Refs. [41, 42]. In 
general, with the mass asymmetry increasing, the max- 
imal dynamical deformation becomes smaller and when 
n approaches to +1, which corresponds to the formation 
of a CN, Bmax approaches to zero. The local excita- 
tion energy becomes larger when the dynamical defor- 
mation develops. For those DNS’s with larger values of 


max 


bBmax, the gain of local excitation energy, 6VpRe(A1) = 
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FIG. 3. (Color online) Maximal values of experimentally 
available evaporation residue cross sections for **Ca induced 
reactions leading to SHN with Z = 112-118 compared with 
the theoretical values calculated at the incident energy in the 
center of mass frame Em. (in MeV) which are indicated in 
the plot. The calculated maximal evaporation residue cross 
sections for °°Ti induced reactions leading to SHN with Z = 
119 and 120 are also shown. For each superheavy nucleus, 
the charge and mass numbers are given and the experimental 
values are shown by black solid squares with error bars and 
the theoretical ones by red open circles. 


Vpes(Aı,t = 0) — Vers(A1,t = oo), is also larger. 
Note that, since Bmax for 7 = +1 is zero, the local 
excitation energy of the CN Héy is always fixed and 
Eéy = Q + Eom. with Q the reaction energy. In Fig. 1 
we plot the potential energy Vpgg (the red dotted curve) 
as a function of 68 for °Ti + 749Cf. It can be seen that 
the potential energy decreases almost linearly with the 
dynamical deformation increasing. Therefore, in order 
to reduce the numerical time, we assume 


SBE) 
Ô Bmax 


Note that a dynamical PES has been calculated micro- 
scopically to describe a continuous transition from the 
initial diabatic potential to the asymptotic adiabatic one 
due to a residual two-body collision [43, 44]. It would be 
an interesting topic to explore the connection between 
the DyPES in the present work and that proposed in 
Refs. [43, 44]. 

With the DNS-DyPES model, we studied systemati- 
cally hot fusion reactions producing superheavy nuclei 
with the charge number Z = 112-120. In Fig. 3, maximal 
values of experimentally available evaporation residue 
cross sections for 48Ca induced reactions leading to SHN 
with Z = 112-118 [4, 45-47] are compared with the the- 
oretical values calculated at the incident energy in the 
center of mass frame Ee.m. which are indicated in the 
plot. The calculated maximal evaporation residue cross 
sections for °°Ti induced reactions leading to SHN with 
Z = 119 and 120 are also shown. For reactions leading to 
SHN with Z = 112-118, the projectile is Ca and tar- 
gets are 8U, 24°Np, 2 Pu, 7 Am, 7° Cm, “Bk, and 
249Cf, respectively. For SHN with Z = 113, 115, and 118, 


Vprs(A1,t) = Vers(A1,t = 0)— 


ôVpps (A1). (10) 
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FIG. 4. (Color online) Evaporation residue cross sections CER 
as a function of the incident energy in the center of mass frame 
Ecm. for the reaction °°Ti+?*°Bk. 


the maximal opp is found in the 3n evaporation residue 
channel and for those with Z = 114, 116, and 117, the 4n 
channel is more favorable. For the element 112, according 
to our calculation, the ER cross section in the 4n chan- 
nel leading to ?8?Cn is a bit larger than that in the 3n 
channel. But the maximal sppr is found in the 3n channel 
in the experiment and in Fig. 3 the calculated and ex- 
perimental results for ?°Cn are given. From Fig. 3, one 
finds a good agreement between the calculation and the 
experiment. We note that the inclusion of the DyPES 
in this work reduces by about an order of magnitude the 
fusion probability and the ER cross section. 

From Z = 112 to 116, the experimental cross section 
staggers and the values are between 1 to 10 pb. For 
Z > 116, the ER cross section decreases almost expo- 
nentially with the charge number increasing. A similar 
trend is also found in the calculated results for SHN with 
Z = 112-118. Whether this trend continues or not is 
a very interesting question. The maximal ER cross sec- 
tions for the SHN 796119 in 5°Ti + 249Bk, 296120 in OTi 
+ *49Cf (the lower circle), and 799120 in °Ti + 751Cf 
(the upper circle) are also shown in Fig. 3. It can be 
seen that, for the element 119, the decreasing tendency 
in the ER cross section continues after Z = 118 and the 
ER cross section is about 0.11 pb. This tendency contin- 
ues also for the element 120 if 74°Cf is used as the target 
and the ER cross section is only about 0.05 pb. How- 
ever, for the 3n-emission channel of °°Ti + 7°4Cf, which 
produces ?98120, the maximal ER cross section is about 
0.25 pb. 

The synthesis of isotopes of the element 119 with °°Ti 
as the projectile has been investigated extensively from 
the theoretical side. As we mentioned earlier, the maxi- 
mal opr from different models or with different param- 
eters vary very much. Let’s take the projectile-target 
combination °°Ti+?4°Bk as an example. The maximal 
opr varies from 35 [17], 50 [16], and up to 570 fb [15]. In 
this work, the production cross sections for the synthe- 
sis of the element 119 are studied with the DNS-DyPES 
model. The excitation functions for °°Ti + ?49Bk lead- 
ing to 794-296119 are represented by solid, dash dot and 
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FIG. 5. (Color online) Evaporation residue cross sections CER 
as a function of the incident energy in the center of mass frame 
Eom. for (a) OTit Cf (b) and OTit" Cf. 


dotted curves in Fig. 4. The maximal cross sections for 
3n and 4n channels are found to be about 0.06 and 0.11 
pb, respectively. 

In an attempt to synthesize SHN with Z = 120 us- 
ing the projectile-target combination 58Fe + ?44Pu, no 
decay chains consistent with fusion-evaporation reaction 
products were observed [5]. According to the sensitivity 
of this experiment, the null result sets an upper limit of 
0.4 pb for ogr. Some predictions have been made for 
the ER cross section for this reaction. For example, the 
maximal opp is predicted to be about 0.1 pb in the 3n 
channel in Ref. [8] and about 0.003 pb in the 4n channel 
in Ref. [17]. In this work, the ER cross section in the 
4n channel is larger and the maximal value is only about 
4 fb which is far below the current experimental limit. 
By examining the excitation function for a more asym- 
metric projectile-target combination, °*Cr + ?4°Cm, we 
find that the maximal cross section is also very small and 
is only about 6 fb in the 4n channel which is similar to 
the results of Ref. [17]. 

Next we investigate °°Ti induced reactions for the syn- 
thesis of isotopes of the SHE with Z = 120. The exci- 
tation functions for °°Ti + 749Cf and °°Ti + 7°'Cf are 
shown in Fig. 5. In both cases the 3n-emission channel 
gives larger ER cross sections than does the 4n channel 
due to the odd-even effects in the survival probability. It 
is found that the maximal cross section for °°Ti + 74°Cf 
is about 0.05 pb. However, for °°Ti + 7°! Cf, the ER cross 
section can be as large as 0.25 pb which is close to the cur- 
rent experiment limit. The reason for a large ER cross 
section in these two reactions is that the fusion proba- 
bility increases considerably with the mass asymmetry 
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increasing; the fusion probabilities for °°Ti + 249:?5Cf 
are about one order of magnitude larger than that for 
58Fe + 744Pu [19]. 

In summary, we developed a di-nuclear system 
(DNS) model with a dynamical potential energy surface 
(DyPES)—the DNS-DyPES model. In this model, the 
development of dynamical deformations of the two nu- 
clei in a DNS is approximately taken into account. This 
is crucial for the determination of the local excitation 
energy of DNS’s involved in the fusion process. With 
the DNS-DyPES model, heavy ion fusion reactions with 
trans-uranium nuclei as targets are investigated. The 
calculated evaporation residue (ER) cross sections are in 
good agreement with available experimental values for 
the reactions producing superheavy nuclei with Z =112— 
118. The projectile-target combination °°Ti+?4°Bk for 
synthesizing the element 119 is studied and the maximal 
ER cross section is found to be about 0.11 pb for the 
A4n-emission channel. For projectile-target combinations 
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